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Abstract. The disruption of a star by the tidal field of a massive black hole is the final outcome of a chain of 
complex dynamical processes in the host galaxy. I introduce the "loss cone problem", and describe the many 
theoretical and numerical challenges on the path of solving it. I review various dynamical channels by which 
stars can be supplied to a massive black hole, and the relevant dynamical relaxation / randomization mechanisms. 
I briefly mention some "exotic" tidal disruption scenarios, and conclude by discussing new dynamical results 
that are changing our understanding of dynamics near a massive black hole, and may well be relevant for tidal 
disruption dynamics. 



1 Outline 

This is a brief, informal review of the stellar dynamical 
mechanisms that play a role in the tidal destruction (TD) 
of stars by massive black holes (MBHs) in galactic nuclei. 
The key issue is to determine how, and at what rate, are 
stars deflected to the extremely eccentric orbits that bring 
them sufficiently close to the MBH for TD. This is the so 
called "loss-cone (LC) replenishment problem", which is 
introduced in its most basic form (single stars around a 
single MBH, scattered singly into the LC) in Sec. [2] to- 
gether with a discussion of the more general framework 
of infall / inspiral processes, and other channels of LC re- 
plenishment. Orbital evolution (denoted here generally by 
the term "relaxation"), whether stochastic or coherent, can 
take various forms and proceed at a wide range of rates 
on different spatial scales. Some lesser known forms of re- 
laxation are presented in Sec. [3] and their impact on TD 
is discussed. TD rate estimates, and the rates of associated 
"exotic" processes, are reviewed in Sec. [4] Sec. [5] explores 
some interesting new dynamical processes that may be rel- 
evant for stars that are formed or captured very close to 
the MBH, and deflected to the LC from tight orbits. Sec. [6] 
briefly summarizes the main conclusions and offers a few 
final comments. 



2 Getting stars to the MBH: the loss cone 
2.1 Basic results 

TD occurs when a star approaches an MBH of mass M. 
closer than the TD radius r, - tf*(M./M*) I/3 , where R+ 
and M+ are the star's radius and mass. Alternative phys- 
ical formulations of the TD criterion are that the star is 
destroyed when its typical density, p* ~ falls be- 

low the density the MBH would have had if its mass were 
spread over the volume r 3 , or that the star is destroyed 
when the crossing time through the disruption zone (r 3 /GM.)' 
falls below the star's free-fall time (Rl /GM*) 1/2 . TD by a 
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MBH can occur only when r, lies outside the event horizon 
at r e = xGM./c 2 , where x = 2 for a non-spinning MBH, 
and can be as small as x — 1 for a maximally spinning 
MBH and a co-orbiting star. Since r,/r e oc p~ ^ 3 'x~ l M~ 2 ^ ', 
it follows that the less massive the MBH, the more com- 
pact the stars that undergo TD can be; for a star with given 
R+, M+, there exists a maximal MBH mass, maxM,, where 
TD is still possible; and the higher the MBH spin, the larger 
is maxM.. 

A star on an orbit with specific energy E = GM./2a 
and specific angular momentum J — J c Vl - e 2 (The Ke- 
plerian limit is assumed, with the convention that E - 
-Zitrue > for a bound orbit, a is the semi-major axis 
(sma), e the eccentricity, and J c = yjGM.a is the specific 
circular angular momentum) will reach periapse inside the 
TD radius, r p = a(l — e) < r t if its angular momentum is 
less than the LC angular momentum J\ c - ylGM.r,(\ + e) 
(Fig. [T]^). All stars that are initially on LC orbits will be 
destroyed in less than an orbital period P. From that time 
on, the rate of TD events will depend on the rate at which 
dynamical mechanisms repopulate these orbits. 

The discreteness of a stellar system guarantees a min- 
imal orbital randomization rate by incoherent 2-body re- 
laxation, on the timescale tz ~ Q 2 P/N i , log Q, where Q = 
M./M+ and N+ia) is the number of stars with sma < a 
(here and below I informally interchange r <-> a). This is 
the time it takes the orbital energy to change by order unity. 
Because the 2-body scattering is local and isotropic (i.e. it 
does not "care" where the MBH is), the time for an order 
unity change in angular momentum can be much shorter 
than t£ if the orbit already has low angular momentum, 
tj ~ [J/J c (E)] 2 tE- It then follows that stars get deflected to 
LC orbits primarily by relaxation of angular momentum, 
and not energy (Fig.[Tp). 

More detailed calculations [e.g. 1 1 indicate that the to- 
tal TD rate is /jd ~ Af*(< r/,)/ (log(J c (a)/Ji c )t E {a)), where 
rh is the MBH's radius of influence, r/, ~ GM./cr 2 , where 
12 cr is the velocity dispersion far from the MBH, and where 
the average (■ • •) is over the volume within r/,. Note that 
the dependence of the rate on the relative size of the LC is 
only logarithmic, so that the TD rate scales approximately 
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Fig. 1. Schematic depiction of LC dynamics in real space and in phase space, and of the difference between infall and inspiral. Left: A 
test star's velocity vector is deflected into the LC by an interaction with a field star (scatterer). When the test star is relatively close to the 
MBH, the angular opening of the LC is large, while the magnitude of a typical accumulated deflection per orbit is small; stars diffuse into 
the LC at a rate much slower than the time it takes a star to reach the MBH, and the LC is almost empty of stars. Conversely, for stars far 
from the MBH, the angular size is smaller than the typical accumulated deflection per orbit, stars are rapidly "kicked" into the middle of 
the LC, which remains almost full. In typical galactic nuclei, most tidally disrupted stars originate from the transition range between the 
empty and full LC regimes, at r ~ r h , where N+M+ ~ M.. Center: The phase space (log/, log E) of infall (plunge) events, such as tidal 
disruption, which happen promptly once the star gets close enough to the MBH. Since tj ~ (J/J c ) 2 t E , once a star reaches a somewhat 
eccentric orbit, there is a substantial probability that it will be scattered into the LC and promptly destroyed. Right: The behavior in 
the presence of a dissipative process, for example the emission of gravitational waves (GW), whose luminosity sharply increases with 
proximity to the MBH. Close enough to the MBH, above some critical energy E crii , (equivalently, inside some critical radius r crit ), where 
the dissipation per period is strong, orbital decay is faster that the 2-body relaxation, and almost all stars eventually inspiral gradually 
into the MBH. Conversely, outside that region, almost all stars ultimately plunge into the MBH. 



as the number of stars in the relevant reservoir, over the 
relaxation time on that scale. 



rates and their parameter dependences. These issues and 
open questions broadly fall in three categories: 



2.2 Plunge vs Inspiral 

TD is one example of an infall process, where the event of 
interest occurs promptly once the star reaches close enough 
to the MBH. It is instructive to consider also the related 
class of inspiral processes, where energy is dissipated grad- 
ually on successive peripassages, and the orbit decays until 
the star is destroyed. The dissipation can be, for example, 



by tidal heating (Sec. 4.3 1, or by the emission of gravita- 
tional waves (GW). In both cases the dissipated power is a 
strong function of the distance from the MBH, and there- 
fore almost all of it occurs near periapse. Since inspiral 
is a race between deterministic dissipation and stochastic 
scattering, there exists a boundary, which is approximately 
constant in energy (equivalently, a typical radius from the 
MBH), that separates infall from inspiral (Fig{TJl). Far from 
the MBH, where the orbital time is long, there is enough 
time between peripassages to be scattered off the nearly 
radial orbit required for a close encounter with the MBH. 
Therefore stars beyond some critical radius r cr j t ultimately 
fall directly (plunge) into the MBH. Conversely, stars in- 
side r cr i t inspiral gradually into the MBH, at a rate /"i nsp ii-ai ~ 
N* [< r cr n(t E )] /(log(7 c /7 fc )f £ ) . This expression, whose in- 
terpretation is similar to that of the infall rate, is decep- 
tively simple, since it is non-linear because r cr ; t depends 
itself on f £ . For GW inspiral in the GC, r cl -i t ~ (9(0.01 pc) 
(Fig. |2jR.). Because typically r cl -i t <sc r/„ and the number of 
stars enclosed within r cr i t is much smaller than in r/,, it is 
generally the case that r inS p ira i ~ <9(0.01)r p i un g e . 



Analytic representation How to achieve a tractable self- 
consistent description of the stellar system? How to set 
boundary conditions correctly? How to describe the dif- 
fusion of orbital elements in phase space? 



Numeric simulation How to simulate the -> oo limit? 
How to deal with the r — > limit (very close encounters; 
General Relativistic (GR) effects)? 



Galactic models Which are the dominant components 
(e.g. one or two MBHs; gas; disks)? What stellar popula- 
tion to assume (mass spectrum)? What stellar density to as- 
sume (distribution function)? What symmetries to assume 
(e.g. spherical, axial)? 

While dealing with questions of analytic representation 
and numeric simulation may require sophisticated analysis 
and state of the art computational techniques, these issues 
are in principle solvable, even if only by brute force com- 
putation. In contrast, the situation with the astrophysical 
uncertainties of galactic modeling presents an irreducible 
problem, since it is unlikely that these could be determined 
with any degree of certainty in the foreseeable future. At 
this time there is no recourse but to explore a large param- 
eter space of plausible scenarios, and map out the range of 
possible outcomes. 



2.3 The challenges 

Many conceptual, technical, and astrophysical challenges 
stand in the way of reliable theoretical predictions of TD 



2.4 Supply channels 

Although much of the theoretical effort of modeling TD 
events focused on the simplest configuration: stars that are 
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singly deflected toward a sing le MBlQ (Sec. fu), there 
are quite a few additional possible channels leading to TD 
events. These include binary stars separated by the tidal 
field of an MBH, or interactions of stars with a circum- 
nuclear gas disk or star disk around a MBH. All of the 
above can also be relevant when the MBH is not station- 
ary in the center of the galaxy, but is recoiling following 
an asymmetric coalescence event, or in the case where the 



stars surround a binary MBH (Sec. 4.2 1. Many of these 
supply channels are discussed in detail elsewhere in these 
proceedings. 



of the Galaxy as a hyper- velocity star 0. A comparison of 
the numbers of bound B-dwarfs observed in the inner 1 
of the GC (the "S-stars", see Fig.|2j^), and the numbers of 
hyper-velocity stars (also mostly B-dwarfs), supports this 
scenario for S-star capture by MP scattering [4|. While the 
capture process itself is "dark", the captured stars are so 
close to the MBH, that a substantial fraction of them are 
likely to eventually undergo luminous TD ||6l|7) (Sec.[5]>. 

3.2 Resonant relaxation 



3 Rapid relaxation mechanisms 

TD rates are determined by the supply channel, and by the 
replenishment rate of the LC. 2-body relaxation is unavoid- 
able, and thus sets a guaranteed minimal replenishment 
rate. However, other processes can contribute, or even dom- 
inate relaxation. Much of the motivation for the study of 
these processes, and most of the empirical tests of these 
ideas, are due to the unique observations of stars in the 
Galactic Center (GC), very close to the Galactic MBH 
(SgrA*) E]|3] (Fig.[2j^). I briefly describe here rapid relax- 
ation by massive perturbers (MPs), by resonant relaxation 
(RR), and orbital evolution in non-spherical systems. 



3.1 Massive perturbers 

the rate of relaxation of a star of mass M* by 2-body en- 
counters with perturbers of number density n p and mass 
M p can be estimated by a simple 'T ~ nvE" collision 
rate calculation. Approximating the relative velocity v by 
the velocity dispersion <x, and the interaction cross-section 
S - jiR 2 with the collision radius R ~ GM p /cr 2 (the clos- 
est approach needed for an 0(1 rad) angular deflection) 
yields a relaxation rate of r oc tipM^/o 3 (integration over 
collisions at larger separations increases the rate by the 
Coulomb factor). The M p dependence of the rate implies 
that even a small density of very massive objects can dom- 
inate 2-body relaxation by stars, when npM 2 ln+M 2 » 1 
(3). This is indeed the case in many galactic nuclei, where 
giant molecular clouds (GMCs), stellar clusters (and pos- 
sibly IMBHs) lower the estimated 2-body relaxation time 
by orders of magnitude relative to that due to stars alone. 

The acceleration of the relaxation rate in the GC, esti- 
mated from observations of GMCs, is dramatic outside r/,, 
where GMCs are found, and substantial also inside r/, (Fig. 
[2p). Since on the r ~ r/, scale, where most of the TD flux 
originates (Sec. [2] Fig. [TJ^), the LC can already be filled 
by stellar 2-body relaxation alone, MPs raise the TD rate 
only by x3 over that expected by stars alone (the effect is 
mostly due to the less massive gas clumps at ~ r/,. This 
conclusion could overturn if IMBHs do exist inside r/,). 

The most dramatic effect of MPs is on the scattering of 
binaries from outside r/, to LC orbits that lead to their tidal 
separation (not to be confused with star-destroying TD), 
resulting in the capture of one of the stars on a tight eccen- 
tric orbit around the MBH, and the ejection of the other out 

1 This dynamical configuration also applies to stars and white 
dwarfs disrupted by intermediate mass black holes (IMBHs), or 
comets disrupted by a central object. 



RR is a rapid process of angular momentum relaxation 
iflOl 13, which takes place in near-spherically symmetric 
potentials where orbital evolution is limited (e.g. fixed el- 
lipses in a Kepler potential; fixed orbital planes in a gen- 
eral spherical potential). In that case, on timescales much 
longer than P, but still shorter than f CO h, the time it takes 
for perturbations to grow and break the coherence of the 
symmetry (e.g. deviations from a Kepler potential due to 
the mass of the stars themselves, or due to GR precession), 
the averaged orbits of the field stars can be viewed as fixed 
mass wires (in a Kepler potential), or fixed annuli (time- 
averaged rosettes in a spherical potential). 

The potential of this near-stationary background mass 
distribution conserves orbital energy, but the finite num- 
ber of stars means that there are order V^* fluctuations 
from isotropy, which give rise to near-constant residual 
torques that rapidly (oc t) change the angular momentum of 
a test star. The accumulated change over a coherence time, 
(AJ)m\\ ~ ( V^V*GM*//?)f co h, where R is the typical size 
of the test star orbit, then becomes the mean free step size 
of an accelerated random walk, which can be written as 
(AJ/J c )(t) = -\Jt/t RR (for t > f co h), where the RR timescale 
is t RR ~ Q 1 P 1 l(N ir t cob ). 

The evolution of the angular momentum then depends 
on the symmetries of the potential and the process that sets 
the shortest coherence time. In a near-Kepler potential, all 
components of the angular momentum are torqued, and 
in particular RR can change the magnitude J, and drive 
a wide orbit into a near radial one that lies inside the LC. 
The relevant coherence-limiting (quenching) mechanisms 
are precession due to the enclosed stellar mass far from 
the MBH, and GR in-plane precession near it (Fig. |2jl). 
This variant of RR is sometimes called "scalar RR". In ad- 
dition, any spherical potential supports a restricted form 
of RR, where the time-averaged mass annuli change only 
the inclination of the orbital plane of a test star, J = J /J, 
but not the magnitude J. In the absence of axial symme- 
try breaking, this is an extremely fast process, since it is 
only quenched by the orbital randomization due to RR it- 
self] The self-quenching timescale f^ h is then set by the 
requirement (/IJ) co h ~ J c , which is a long timescale, f^ h ~ 
QP/ "vW*- This variant of RR is sometimes called "vector 
RR". 

RR is easily detected in A^-body simulations ifTTIl . and 
there is also some evidence that it is operating in the GC [9 1 
(Fig. |2jl). Its direct effect on the integrated TD rates is not 
large (up to x2 enhancement [10]), since most of the TD 
flux is from r ~ r/,, where Af+M* ~ M. and orbital preces- 
sion due the enclosed stellar mass quenches RR. However, 

2 The separation of background stars / test star is arbitrary-all 
the stars in the system are simultaneously torqued by RR. 
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Fig. 2. Lessons from the GC. Left: A not-to-scale schematic of the various dynamical components that are observed in the GC around 
SgrA*, or are hypothesized to exist there (marked by [. . .]) |8|. Note in particular the massive GMCs and the circum-nuclear molecular 
disk at or outside the MBH's radius of influence; the mixed old and young stellar population inside r h (only the giants and B-dwarfs 
are above the present-day detection threshold); the (9(100) blue giants that are concentrated in one (or perhaps two) coherently-rotating 
warped disks, which extend inward to a sharp inner cut-off at ~ 0.04 pc (~ 1"); and the spectroscopically and dynamically distinct 
population of ~ 40 B-dwarfs on isotropic orbits inside ~ 1" (the "S-cluster"). Center: The GMCs and molecular clumps of the circum- 
nuclear disk dramatically shorten the relaxation time outside r h (blue and green lines, for two limiting assumptions about the masses of 
the GMCs), relative to that by stars alone (red line), and also reduce from afar the relaxation time inside r h (circles) |4 8 1. Right: Possible 
signature of RR in the inner pc of the GC [9]. Significant relaxation is expected where the relaxation times are shorter than the typical 
ages of the different dynamical component. 2-body relaxation is faster than scalar (/-changing) RR only far from the MBH, where the 
enclosed stellar mass becomes substantial, or extremely close to the MBH, where GR precession is strong. The two lines for scalar RR 
represent different assumptions about the mean stellar mass. Scalar RR is an important process for the dynamics of GW sources, and 
may be also implicated in the randomization of eccentricities of the S-stars (6) [7) and the relaxed giants. Vector (/-changing) RR is 
everywhere faster than 2-body relaxation, as long as the average potential is close to spherical, Vector RR may explain the inner cutoff 
of the young star disk, and the randomization of the orbital inclinations of the relaxed giants, (reproduced with permission from the 
Astrophysical Journal). 



RR may play a more significant role for TD events by stars 
on tight orbits (Sec. [5]). 



3.3 Non-spherical systems 

A high level of symmetry in the stellar potential is associ- 
ated with conserved quantities, which restrict orbital evo- 
lution and can slow down LC refilling to the minimal rate 
provided by 2-body relaxation (a collisional effect outside 
the domain of potential theory). Conversely, less symmetry 
leaves more freedom for orbits to evolve, and opens new 
possibilities for efficient replenishment of loss-cone orbits. 

Axisymmetric systems have centrophilic (center appro- 
aching) orbits, which lead to a x2 enhancement in the TD 
rate over that predicted for spherical systems [12|. An- 
other consequence is the higher incidence of deep TDs, 
with very small periapse distance, which may result in tidal 



detonation (Sec. 4.3 i. Triaxial systems can support a large 
fraction pf chaotic orbits, which can lead to a dramatic 
xlO - 100 enhancement while the chaotic orbits persist 
Ifl3l . This conclusion depends however on whether chaotic 
orbits can exist in the presence of the stabilizing influence 
of the central MBH. 



4 Tidal disruption rates 
4.1 Single MBH in steady state 

The basic configuration to consider is TD on a single MBH 
surrounded by stellar cluster in steady state. TD rate esti- 
mates for such a configuration have a long history, starting 
in the mid 1970's. Early analyses focused on the similar 
process of a star disrupted by an IMBH in a globular cluster 
lfT4l[T5ll . The first estimate for TDs on a MBH in a galactic 
nucleus yielded a very high rate , of 10~ 3 - 10 M Q yr~' gal -1 



fl6l . Contemporary rate estimates, with more careful mod- 
eling of the process, yield consistently lower values. These 
are summarized in table Q] 



4.2 Binary MBHs beyond steady state 



The near ubiquity of nuclear MBHs from very early cosmic 
times, and the bottom-up structure formation paradigm, to- 
gether imply that MBH pairs and bound MBH binaries 
must be a commonly occurring, probably transient phase 
in galaxy evolution. Recent studies have found that hard, 
high mass ratio binary MBHs are efficient dismptors. This 
can come about by Kozai forcing of stellar orbits around 
the primary (more massive) MBH, induced by the sec- 
ondary MBH [24 1, or by chaotic mixing due to the sec- 
ondary |25l . These effects can lead to a transient phase 
of very high TD rates, < 1 yr~', and account overall for 
(9(10%) of the steady state cosmic TD rate. Another con- 
sequence of this high TD rate is the possibility of multiple 
events from the same galaxy over an observationally rele- 
vant short timescale, which could provide a signature of a 
binary MBH (26). 

Certain combinations of orbital angular momentum and 
MBH spin orientations lead to a strong asymmetry in the 
emitted GW flux and a very fast recoil velocity for the 
newly coalesced MBH [27|. The recoiling MBH retains 
a small cluster of tightly bound stars with it. The recoil 
can promptly refill the LC and briefly raise the TD rate to 
(9(0.1 yr~'), thereby providing a prompt EM counterpart to 
the GW signal of the coalescence [28]. Much later, as the 
recoiling MBH moves farther away from the center of the 
galaxy, the small dense cluster it carries with it can provide 
off-center or even extra-galactic TD events at rates between 
0.01 [29) to < 1 El of those of non-recoiling MBHs. 
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Table 1. Contemporary TD rate estimates for a single MBH in steady state 
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4.3 TD-related exotica 

It is interesting to note that there are exotic processes that 
must occur, and at substantial rates, in association with TD. 
These happen when the orbital periapse is either just out- 
side the tidal disruption radius r„ or very deep inside it. 



Near misses Stars that are deflected from a tight enough 
orbit to an eccentric one with periapse within ~ 2r, can en- 
ter an inspiral phase, similar to that of GW inspiral events, 
as "squeezars" [31], stars that brightly shine by the orbital 
energy dissipated as they tidally interact with the MBH. 
This inspiral-type mechanism has a rate that is typi cally ~ 
0.05 of the TD rate, 0.1-1 yr 1 in the GC (Sec.|Z2|. When 
stars are not captured on an inspiral orbit, they undergo a 
single very strong tidal scattering interaction 021 . which 
can spin them up and mix them, thereby affecting their sub- 
sequent spectroscopic signature and evolution. Such event 
are predicted to occur at 0(1) of TD rate, 10 4 - 10 s yr in 
GC. Finally, the extended, low-gravity envelopes of giant 
stars are susceptible to tidal stripping, leaving behind the 
hot bare cores [33 1. These events are estimated to occur at 
~ 0.01 of the TD rate in the GC. 



Direct hits When the periapse of the incoming star is 
<k r,, and the stellar core is in a state that is favorable 
to fast nuclear burning processes (post-hydrogen burning), 
the strong tidal pancake compression the star experiences 
as it crosses periapse can lead to a supernova-like tidal det- 
onation event [34, 35] . The rate of such events is estimated 
at < 0.1 of the TD rate. 



5 Tidal disruptions from tight orbits 

Standard incoherent 2-body relaxation results in loss-cone 
refilling that is dominated by stars from typical radii r ~ r/,, 
where the dynamical timescale is long, the MBH poten- 
tial no longer dominates (and therefore RR is inefficient), 
and GR effects are completely negligible. However, there 
are compelling reasons to believe that there are situations 
where stars on much tighter orbits can contribute signifi- 
cantly to the total TD rate. This can be the case, for ex- 
ample, when there is substantial in-situ star formation in 
gravitationally unstable accretion disks, as is observed in 



the GC on the (9(0. lr/,) scale [36|; or for stars captured 
around the MBH by the tidal separation of binaries |37l 
(Sec. [XT); or for stars that are tightly bound to a recoiling 

mbhooIED. 

In such cases, dynamical mechanisms that have not 
been considered before may come into play. I focus here on 
two newly uncovered effects: stellar relaxation by anoma- 
lous diffusion [38|, and the Schwarzschild Barrier (SB) in 
phase space l39l . While it is too early to tell whether these 
effects have a significant impact on TD events, their recent 
discovery highlights the possibility that basic key compo- 
nents may still be missing from our understanding of dy- 
namics near MBHs. 



Anomalous diffusion Normal diffusion can be viewed 
as a random walk process, where the accumulated change 
AE or A J grows with time as oc yft (I omit here for brevity 
discussion of drift). This is a consequence the Central Limit 
Theorem (CLT) applied to multiple scattering events. The 
CLT is valid since each of those is physically limited to 
some finite range of AE or A J, and is thus well-defined 
statistically. However, the rate of convergence to the CLT 
can be very slow, depending on how nearly-divergent is 
the physics of the scattering. In that case, on timescales 
shorter than CLT convergence, the evolution of E or / is 
by anomalous diffusion. The rate can be very different from 
\t; it may depend on AE itself (non-self-similar evolution); 
it can result in improbably high rates of "high-cr" events 
(when mis-interpreted as Normal probability events). The 
Newtonian gravitational force with its formal divergence 
at r\2 — > is exactly such a special case (Fig. [3p). On 
timescales shorter than ~ 0.1 f£ there are substantial devia- 
tions from the expected statistics of normal diffusion. 

Anomalous diffusion may be particularly relevant close 
to the MBH, where it could change the branching ratios 
between infall and inspiral (effectively changing r cr jt, see 



Sec. 2.2 1. This can be seen by noting that inspiral can oc- 
cur only when the dissipation timescale is shorter than the 
relaxation time, f^iss ~ {E/AEa^P < tj(J/ c ), where AE& %% 
is the energy dissipated per peripassage, and where both 
zfEdiss and tj are evaluated close to the LC boundary. 

Since t^ ss /tj(Ji c ) oc A^* oc a 3 ~ a for an n* oc r~ a cusp 
(Sec. |2.1[ ), then if anomalous diffusion in A J (not yet ex- 
plored) leads to overall faster /-relaxation, it will shift r cr ; t 
to a smaller radius, implying more prompt TDs and fewer 
tidal inspirals than estimated for normal diffusion. 
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Fig. 3. Previously unconsidered dynamical effects that could modify TD rates from orbits near the MBH. Left: Anomalous diffusion on 
timescales shorter than t E leads to faster than oc yi evolution in phase space (crosses from yV-body simulation and theory in red trace the 
relative change in energy AE/E), and in particular makes "high-cr" events much more common than they would be in normal diffusion 
(Gaussian probability distribution function in green) |38|. Center: The phase-density of 50 M Q stellar mass "black holes" around a 
M. = 10 6 M G MBH in (log[l - e],loga/l mpc) ~ ( log y, log £) phase space, as obtained by Monte Carlo experiments [ 40 1. In the 
absence of GR, RR rapidly throws all stars into the MBH across the last stable orbit (LSO). Right: With GR in-plane precession, the 
Schwarzschild Barrier (a resonance between the RR torques and the GR precession) 1 39 1 deflects stars back to lower eccentricity orbits 
and allows gradual GW inspiral to take place (see tracks with nearly constant r p parallel to the LSO line). 



The Schwarzschild Barrier RR torquing of orbits (Sec. 
|3.2| Fig |2jl) becomes more important closer to the MBH, 
until it is quenched by GR in-plane precession very close 
to it. It was hypothesized that this quenching is crucial for 
"saving" stars from plunging into the MBH, and allowing 
them to inspiral gradually (Fig. |3p). Relativistic A^-body 
simulations ll39ll revealed that the GR precession not only 
quenches RR and enables inspiral, but in fact creates a re- 
flective barrier in phase space, which can be approximately 
understood as arising from a resonance between the GR 
precession frequency and the rate of angular momentum 
change due to the RR torques (Fig.[3jl). Stars can cross the 
SB only by 2-body scattering. Once past the barrier, they 
are free to inspiral. Current analysis of the SB phenomenon 
has focused only of GW inspiral. The implications for TD 
events from tight orbits are still unknown. 

6 Summary 

It is clear that there are multiple plausible channels by 
which stars can be deflected into LC orbits that lead to 
TD events. However, quantitative estimate of the cosmic 
rates of such events are severely limited by astrophysical 
uncertainties. On the theoretical side, there is a clear need 
for a better understanding of the conventional dynamical 
mechanisms that are known to play part in the replenish- 
ment of the LC, and there are also new, recently discovered 
mechanisms that need to be taken into account. In spite of 
these uncertainties, it seems quite realistic to expect that 
TD rates should be high enough to be observationally in- 
teresting. Indeed there are multiple efforts to approach the 
problem empirically, and detect TD events (see in these 
proceedings). 

As is evident from this review, the dynamics of TD 
events are tightly linked to those of GW sources through 
the infall / inspiral connection. The study of one bears di- 
rectly on the other. One additional important connection is 
due to the unfortunate fact that the electromagnetic signa- 
ture of TD events (discussed elsewhere in these proceed- 
ings) may not be distinct enough from other types of AGN 



variability to offer definitive identification. This may ulti- 
mately be achieved only by the detection of low -frequency 
GW counterparts to TDs. 

This work was supported by ERC Starting Grant No. 202996 and 
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